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The colour fields, created by a static gluon-quark-antiquark system, 
are computed in quenched SU(3) lattice QCD, in a 24 3 x 48 lattice at 
f3 = 6.2 and a = 0.07261(85) fm. We study two geometries, one with a U 
shape and another with an L shape. The particular cases of the two gluon 
glueball and quark- ant iquark are also studied, and the Casimir scaling is 
investigated in a microscopic perspective. This also contributes to under- 
stand confinement with flux tubes and to discriminate between the models 
of fundamental versus adjoint confining strings, analogous to type-II and 
type-I superconductivity. 
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1. Introduction 

In this paper, we present a value for the the dual gluon mass in a SU(3) 
lattice QCD gauge independent and a detailed study of the Casimir scaling. 
In section II, we introduce the lattice QCD formulation. We briefly review 
the Wilson loop for this system and show how we compute the colour fields 
and as well as the lagrangian and energy density distribution. In section 
III, the numerical results are shown. We present results for the colour field 
profiles in the mid flux tube section for the static hybrid gqq, in a U shape 
geometry. A detailed study of the Casimir scaling is done and we present a 
value for the effective dual gluon mass and some values found in literature 
for the effective dual gluon mass and gluon mass. Finally, we present the 
conclusion in section IV. 
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2. The Wilson Loops and Colour Fields 

The Wilson loop for the static hybrid gqq was deducted in [2 EJ HJ [3], 
therefore we only present the fundamental expressions. The Wilson loop for 
this system is given by 



W, 



gqq 



Wl w 2 _ _ Ws 



(1) 



where Wi, W2 and W3 are the simple Wilson loops shown in Fig. [Tj 




(b) 



Figure 1: (a) Wilson loop for the gqq and equivalent position of the static 
antiquark, gluon, and quark. |(b)| Simple Wilson loops that make the gqq 
Wilson loop. 



In order to improve the signal to noise ratio of the Wilson loop, we use 
APE smearing, [3], with w = 0.2 and iterate this procedure 25 times in the 
spatial direction. To achieve better accuracy in the flux tube, we apply the 
hypercubic blocking (HYP) in time direction, [5], with a\ — 0.75, = 0.6 
and as = 0.3. Notice that we only apply the smearing technique to the 
Wilson loop. 

The chromoelectric and chromomagnetic fields are given by, 



{Ef) = (Poi) 



(WP 0i ) 
(W) 



(2) 
(3) 



where the jk indices of the plaquette, P, complement the index i of the 
chromomagnetic field. The energy {%) and lagrangian (£) densities are 
given by 



(4) 
(5) 
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3. Results 

Here we present the results of our simulations with 286 24 3 x 48, f3 = 6.2 
quenched configurations generated with the version 6 of the MILC code 
[6], via a combination of Cabbibo-Mariani and overrelaxed updates. The 
results are presented in lattice spacing units of a, with a = 0.07261(85) fm 
or a" 1 = 2718 ± 32 Me^. 



va y 




(a) U shape geometry. (b) L shape geometry. 
Figure 2: gluon-quark-antiquark geometries, U and L shapes. 

In this work two geometries for the hybrid system, gqq, are investigated: 
a U shape and a L shape geometry, both defined in Fig. |2j In the L shape 
geometry only the case when the gluon and the antiquark are superposed, 
the quark-antiquark case, is studied. The use of the APE (in space) and 
HYP (in time) smearing allows us to have better results for the flux tube, 
Fig. |3j while suppressing the fields near the sources. 
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Figure 3: Comparison between results for lagrangian density with and with- 
out Hypercubic Blocking (HYP) smearing in time. 



3.1. U profiles and Casimir scaling 

In Fig. [4] we present the profiles for the U geometry for / = 8 and d 
between and 16 at y = 4. We can see the stretching and partial splitting 
of the flux tube in the equatorial plane (y = 4) between the quark and 
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antiquark. For d — 2 and 4 at y — 4 the separation between the two flux 
tube is not visible, this is due to the overlap in the tails of the flux tube 
which contributes for the total field, and for large separations the tails of 
the flux tubes contributes to a non zero field at x = 0. 

We measure the quotient between the energy density of the two gluon 
glueball system and of the meson system, in the mediatrix plane between the 
two particles (x = 0). The results are shown in Fig. [5j In Fig. 5a we present 
the results for r = y at x = z = and in Fig. 5b we present the results for 
r = (x, z) at y = 4. We make a constant fit to the data in Fig. [5| the result 
for Fig. [talis 2.25096 ± 0.0244972 and for Fig. [5b] is 2.23591 ± 0.0598732. 



As can be seen, these results are consistent with Casimir scaling, with a 
factor of 9/4 between the energy density in the glueball and in the meson. 
This corresponds to the formation of an adjoint string. 





(a) Chromoelectric Field 



(b) Chrome-magnetic Field 





(c) Lagrangian Density (d) Energy Density 

Figure 4: Results for the U geometry at y = 4 and z = 0. 



3.2. Dual Gluon Mass 

In 1970's, Nambu |7J, 't Hooft |8J and Mandelstam [9] proposed an inter- 
esting idea that quark confinement would be physically interpreted using the 
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(a) r — y at x — and z — (b) r = (x, y) at y = 4 

Figure 5: Casimir scaling. 



dual version of the superconductivity, the QCD vacuum state to behave like 
a magnetic superconductor. The chromoelectric field originated by a qq pair 
is squeezed by Meissner effect into a dual Abrikosov flux tube, giving rise 
to the confining linear potential, the field is confined into flux tubes, QCD 
strings. Color confinement could be understood as the dual Meissner effect. 
In common superconductivity the magnetic field decays with B ~ e ~ r / x L 
and this could be interpreted in terms of an effective mass for the photon 
m 7 = 1/Xl- There also is evidence for the dual superconductor picture from 
numerical simulations of QCD, some studies have point a similar behavior 
in QCD, [mEll 02]. 

We tested two functions, ae" 2/if and clKq (fir), where \i — is the 

penetration length and Ko the modified Bessel function of order zero. So, 
in this case we have \i as the dual gluon mass. Fitting the chromoelectric 
field and the lagrangian density section in the mid distance of the flux tube 
of the meson and the two gluon glueball, we obtain the results presented in 
Table [T] for the effective dual gluon mass, of the order of ~ 1 GeV. Some 
values found in literature, for the effective dual gluon mass, [T3 l H^ f [T5 | fT6| . 
and the effective gluon mass, [TT] , point for a same value. 



4. Conclusions 

When the quark and the anti-quark are superposed, this corresponds to 
the formation of an adjoint string between the two gluon and agrees with 
Casimir Scaling measured by Bali |18| . This can be interpreted with a type- 
II superconductor analogy for the confinement in QCD with repulsion of 
the fundamental strings and with the string tension of the first topological 
excitation of the string (the adjoint string) larger than the double of the 
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ae~ 2 ^ r 


a Kq (/i r) 


/i (GeV) 




/i (GeV) 


X 2 /dof 




1.170 ±0.228 


1.069 


0.805 ±0.287 


1.827 


£(l)(a) ( r ) 


1.170 ±0.119 


0.512 


0.865 ±0.188 


1.203 


^(2)(a) W 


1.231 ±0.286 


1.547 


0.881 ±0.334 


2.084 




1.210 ±0.056 


0.887 


0.897 ±0.085 


1.185 


c ww ( r ) 


1.208 ±0.068 


0.560 


0.909 ±0.099 


0.909 


^(2) (6) ( r ) 


1.210 ±0.063 


1.162 


0.889 ± 0.097 


1.262 


£(2)(6) (r) 


1.191 ±0.031 


1.066 


0.899 ± 0.048 


1.106 



Table 1: Results for the dual gluon mass, where (1) is for the two gluon 
glueball and (2) for the quark-antiquark cases, and (a) at y = 4 and z = 
with r — x and (b) at y = 4 with r — (x, z). 



fundamental string tension. 

We present a value for the dual gluon mass of ^ 1 GeV which is gauge 
independent. 
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